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ABSTRACT 

We present deep WFC3 grism observations of the candidate z ~ 12 galaxy UDFj -39546284 in the HST Ultra 
Deep Field (UDF), by combining spectroscopic data from the 3D-HST and CANDELS surveys. The total 
exposure time is 40.5 ks and the spectrum covers 1.10 < A < 1.65 /im. We search for faint emission lines 
by cross-correlating the 2D G141 spectrum with the observed //i6o morphology, a technique that is unique to 
slitless spectroscopy at HST resolution. We find a 2.1 a detection of an emission line at 1.599 /im — just redward 
of the JHi/^Q filter — with flux 3.5 ± 1 .3 x 10~^^ erg s~^ cm~^. Assuming the line is real, it contributes 1 10 ± 40% 
of the observed Hi^q flux and has an observed equivalent width > 7300 A. If the line is confirmed, it could be 
Ly-a at z = 12.12. However, a more plausible interpretation, given current results, could be a lower redshift 
feature such as [O III]A4959,5007 at z = 2.19. We find two other 3D-HST [O III] emitters within 1000 km s"^ 
of that redshift in the GOODS-South field. Additional support for this interpretation comes from the discovery 
of a bright "[O III] blob" with a secure G141 grism redshift of z = 1.605. This object has a strikingly large 
observed equivalent width of nearly 9000 A that results in similar "dropout" colors as UDFj -39546284. 
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1. INTRODUCTION 

Over the past few years deep imaging programs with the 
WFC3 camera on the Hubble Space Telescope (HST) have 
opened up the Universe at z > 8 for systematic study (e.g., 
Bouwens et al. 2011a, 2012a; Oesch et al. 2012; ElHs et al. 
2013; Coe et al. 2013). Photometric data can provide strong 
constraints on the redshifts and properties of very high red- 
shift galaxies (e.g., Bouwens et al. 2011b, 2012b). How- 
ever, the interpretation is hampered by the lack of spectro- 
scopic information that provides unambiguous redshift deter- 
mination and the contributions of emission lines to the broad 
band fluxes. The equivalent widths of rest- frame optical lines 
increase with redshift (e.g., Fumagalli et al. 2012), and the 
emission line contribution can be significant even in the broad 
IRAC filters (e.g., Labbe et al. 2010; Ono et al. 2010). 

Obtaining high quality spectra of these very distant ob- 
jects is a formidable challenge because of their faintness and 
because key emission and absorption features are redshifted 
into the observer's near-IR. Apart from a gamma ray burst 
at z = 8.1 (Salvaterra et al. 2009; Tanvir et al. 2009) the only 
claimed (and disputed; Bunker et al. 2013) detection is a Ly-a 
line at z = 8.6 in a z-dropout galaxy in the UDF (Lehnert et al. 
2010). 

Here we describe a deep WFC3 G141 grism spectrum of 
one of the most intriguing high redshift candidates yet dis- 
covered: UDFj-39546284. This object was identified as the 
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first robust z ^ 10 candidate by Bouwens et al. (2011a): it 
was only detected in ifigo in the HUDF09 dataset (see also 
Oesch et al. 2012). Recent observations from the HUDF12 
program (Ellis et al. 2013) confirm the ifi6o detection and 
show that it is undetected in ///140 (Ellis et al. 2013; Bouwens 
et al. 2013). The > 1.4 mag dropout (2cr) in the JHuq filter 
pushes its likely redshift to z ~ 12, although alternative expla- 
nations cannot be ruled out (Ellis et al. 2013). Here we report 
a 2.7(7 detection of an emission line for UDFj-39546284 that 
would account for most or all of its /fi6o flux. Though this 
line could be Ly-a at z = 12, we argue that it is more likely 
to be [O III]A4959,5007 at z = 2.19 based in part on the dis- 
covery of a bright low-redshift analog with one of the highest 
[O III] equivalent widths ever observed (8700 A) and broad- 
band colors similar to UDFj-39546284. Magnitudes are in 
the AB system throughout. 

2. OBSERVATIONS & ANALYSIS 

The UDF was observed with the WFC3/G141 grism as part 
of the 3D-HST survey (GO-12177; PI: van Dokkum; Bram- 
mer et al. 2012a). The observations comprise four sepa- 
rate visits at the same orientation with G141 integrations of 
4.7 ks (two orbits) each. The UDF was also observed by 
the CANDELS supernova follow-up program (GO- 12099, PI: 
Riess) when a supernova at z = 1.55, nicknamed "Primo", was 
discovered in the first epoch of CANDELS imaging in the 
GOODS-South field(see Rodney et al. 2012). The field was 
observed in two visits on 2010 Oct 26 and 2010 Nov 01 with 
G141 integrations of 6.6 and 15 ks, respectively. The posi- 
tion angle of the grism dispersion axis is rotated by 6.4 de- 
grees between the two visits, which together are rotated 50 
degrees with respect to the 3D-HST visits. The relation of the 
G141 pointings to the deep HUDF09+HUDF12 imaging area 
is shown in Fig. 1. 

The Primo and 3D-HST visits were pre-processed sep- 
arately and identically using the threedhst pipeline^ 

^ http : //code . google . com/p/threedhst/ 
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(Brammer et al. 2012a). The individual exposures of each 
visit were then combined into a mosaic by interlacing the 
pixels to a subgrid of pixels that are exactly half their orig- 
inal size. Interlacing pixels from input images which have 
0.5 pixel offsets (by design in 3D-HST and rounded for 
CANDELS-Primo) results in a one-to-one correspondence 
between input and output pixels (e.g., van Dokkum et al. 
2000) and offers the key benefit of preserving the individual 
pixel errors. 

With the HUDF12 Hi^q image defining the 1 .6 /im flux den- 
sities and spatial morphologies, we compute a full quantitative 
contamination model of all objects in the field to Hi^q < 28. 
The modeling technique, first presented in Brammer et al. 
(2012b), includes the first 4 spectral orders, among them the 
compact zeroth order that can resemble emission lines. The 
models of bright objects with //i6o < 23 are refined based on 
their observed spectra, while fainter objects are assumed to 
have continua flat in units of fx. Therefore, emission lines in 
fainter galaxies are not included in the contamination model. 

To generate a deep 2D spectrum from all UDF grism inte- 
grations, we first extract two-dimensional spectra from each 
visit. These spectra are aligned in two dimensions to the near- 
est pixel in the subsampled grid and the aligned pixels are then 
combined with inverse variance weights including a term to 
down- weight contaminated pixels. Stacked spectra from the 
central 3 arcmin^ of the UDF have a combined exposure time 
of 40.5 ks. 

3. A TENTATIVE EMISSION LINE IN UDFj-39546284 
3.1. The Deep Near-IR Spectrum 

The two-dimensional 40.5 ks G141 spectrum of UDFj- 
39546284 is shown in Fig. 2. While some contamination 
remains in the stacked spectrum, the orientations of the 
UDF/G141 visits turn out to be nearly perfect for avoiding sig- 
nificant contamination of the UDFj-39546284 spectrum, par- 
ticularly at A > 1.4/im. The contamination-subtracted spec- 
trum is remarkably clean with none of the systematics that 
typically plague ground-based NIR spectroscopy (see, e.g., 
Lehnert et al. 2010). The cluster of bright pixels at the top 
of the spectrum is likely a faint contaminant emission line not 
included in the model; it lies well away from any contribution 
from UDFj-39546284. 

The "salt-and-pepper" appearance of the 2D spectrum is 
the result of the uncorrected noise of the individual pix- 
els. The distribution of cleaned pixel values is gaussian with 
(Jobs = 0.0033 electrons s~^ Random gaussian deviates scaled 
by the propagated pixel errors have cr^^ = 0.0036 electrons s~^ . 
This slight overestimation of the errors by 10% is significantly 
lower than the factors of 60-70% that result from typical ap- 
plications of the Drizzle algorithm (Fruchter et al. 2009). 

3.2. A 2a Emission Line Detection at X= 1.599 /im 

Slitless HST grism spectra offer a distinct advantage com- 
pared to ground-based data for detecting weak emission line 
features: the 2D spectrum provides an image of the galaxy at 
the wavelength of the emission line (see, e.g., Brammer et al. 
2012b), and the morphology of the galaxy is known a priori 
from the deep broad-band WFC3 imaging. As long as the 
emission line morphology and the broad-band morphology 
are similar (which is generally the case at z ~ 1 ; Nelson et al. 
2012), or the line contributes significantly to the broad-band 
flux (which is the case here; see below), we can search for 
features in the 2D spectrum whose morphology corresponds 
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Figure 1. Exposure map of the 3D-HST and CANDELS G141 grism point- 
ings in the HUDF. The vectors indicate the direction and extent of the first- 
order G141 spectra. The locations of the CANDELS supernova and the two 
faint UDF galaxies described below are indicated. These objects lie in the 
area of deepest G141 grism coverage, with a total exposure time of 40.5 ks 
(17 orbits). 

to that of the broad-band image. 

While there are no obvious emission lines in the raw 2D 
spectrum, there is a clump of pixels at ~ 1.6/im which is 
marginally enhanced compared to its surroundings. To im- 
prove the signal-to-noise of this feature, we cross-correlate 
the 2D spectrum with a kernel defined by the central R = 0i^3 
of the deep Hm thumbnail of UDFj-39546284. This kernel 
is constructed by extracting kernels for each of the individual 
UDF/G141 visits with different orientations and combining 
them weighting by the median error in their corresponding 2D 
grism spectra. Thus, the object profile is slightly smoothed but 
will reflect the 2D morphology of lines in the stacked spec- 
trum. The cross-correlation spectrum is shown in the right- 
hand panels of Fig. 2, and the strength of the feature at 1.6 fim 
is enhanced compared to the raw 2D spectrum. We verified 
that the weak feature is not visible in a similar analysis of any 
single UDF/G141 visit, confirming that the feature does not 
arise from a flux excess in one of them individually, such as a 
group of hot pixels or an un-flagged cosmic ray. 

We show a one-dimensional spectrum extracted along the 
trace in Fig. 3. We compute an associated uncertainty at each 
pixel along the trace, highlighted in blue, by cross-correlating 
the squared error array with the same Hieo kernel. Further- 
more, we compute a 2D model spectrum as in Fig. 2 with 
known position and integrated line flux, and extract its cross 
correlation spectrum in the same way as the observed spec- 
trum. The feature at A = 1 .599 ± 0.004 /im is detected at 2.1a 
with an integrated line flux of 3.5 ± 1.3 x 10~^^erg s~^ cm~^. 
The probability of finding a gaussian noise feature with equal 
or greater significance is ~10% for 30 independent resolution 
elements (-3000A). 

3.3. The Reality of the 1.599 fim Feature 

There are regions in the 2D cross-correlation spectrum 
with apparently similar significance to the feature at 1.6 /im. 
Many of them can be rejected because they do not fall pre- 
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Figure 2. The deep G141 grism spectrum of UDFj-39546284. The left panels show the pixels as observed and modeled, and the right panels show the spectra 
cross-correlated with the //i6o object profile to enhance weak features. The first row shows the stacked spectrum as observed. The second row shows the 
contamination model (§2), with some aliasing due to the unequal weights of the interlaced Primo mosaics. The third row shows the observed spectrum cleaned 
of this contamination. There is a weak emission feature in the cross-correlation spectrum at A = 1.599/im. The bottom row shows the emission line model of 
UDFj-39546284 where the morphology is derived from the observed Hi^q image (see also Fig. 3). 





UDFj-39546284 
Lya,z= 12.12? 
[OIII]A.5007,^=2.19? 





F160W kernel ^ Contam. 

Uniform kernel (R=0.30") 




UDF-40106456 
(H|3H[0III]-hHa,^= 1.303 



1.2 1.3 1.4 1.5 

A / jUm 



1.6 



Figure 3. Optimal extraction of the UDFj-39546284 spectrum (left) and that of a comparison object, UDF-40106456(right). The top panels show the stacked 2D 
spectra and the spectra cross-correlated with the Hk^q image kernel, which is indicated at left. The upper spectra show the 2D cross-correlation extracted along the 
trace and scaled by the G141 sensitivity curve and A A of the spectrum, which for a uniform kernel would give an integrated line flux with the units as labeled. The 
same cross-correlation is applied to the 2D pixel variance array, yielding the uncertainties shaded in gray. The bottom panels show the signal-to-noise spectrum. 
While there are additional apparently significant features at A < 1.5 jim, these are likely due to contamination rather than coming from UDFj-39546284 (§3.3). 
The gray curves in the bottom panels show the S/N 'mR = Of^3 apertures along the trace. The red curves indicate the potential contamination from nearby sources. 
The filled orange regions are model emission line spectra, whose shapes match those of the observed lines. 
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cisely along the trace of the grism spectrum. There are en- 
hancements found near the trace at A = 1 . 14 /im and 1 .22 /im, 
though these wavelengths suffer from higher contamination 
perhaps consisting of faint emission lines not included in the 
contamination model. The HSTfW¥C3 broad-band photome- 
try places strong constraints on the possibility that these lines 
come from UDFj -39546284: either of the bluer lines alone 
would result in Fios or Jus two magnitudes brighter than the 
2cr HUDF12 limits. 

In contrast, the 1.599/im line is consistent with the 
HUDF12 photometry. Within the same 0'/3 aperture used to 
measure the line flux, we measure //i6o = 29.04 for the //i6o 
kernel. This is 27% brighter than the Ellis et al. (2013) mag- 
nitude measured within an Of'25 aperture, but is only 73% 
of the total magnitude measured by Bouwens et al. (2013). 
These differences are consistent with aperture corrections of 
the extended source, which shows a faint extended tail to the 
NE in the //i6o image (Bouwens et al. 2013). A pure emission 
line with the observed flux would result in Hieo = 28.92; the 
line accounts for 110 ±40% (la) of the Hi^o flux. This cor- 
responds to a (l,2)cr limit on the observed- frame equivalent 
width of the line of EW > (7300, 1300) A. 

We test our technique by analyzing faint lines of another 
UDF galaxy with a secure redshift, UDF-40106456(right pan- 
els of Fig. 3). This galaxy has //i6o = 27.4 and unambigu- 
ous emission lines of Ua and [O III]A4959,5007 at z = 1.303 
visible in the stacked 2D spectrum before cross -correlation. 
Using the cross-correlation technique we find a significant 
detection of H/3 even though the line is barely visible in 
the original spectrum. Two spurious features are detected at 
1.24/im and 1.29 /im with comparable flux as the tentative 
line in UDFj-39546284. These are both clearly associated 
with residual contamination from another source at the top 
of the 2D spectrum. The Ua line flux of UDF-40106456 is 
9.5 ± 1.0 X 10~^^erg s~^ cm~^, just three times brighter than 
the UDFj-39546284 emission line. The detection significance 
is S/N'^10 and the cross-correlation technique enhances the 
detection significance by a factor of two compared to simple 
photometry along the trace within an equivalent Of' 30 aper- 
ture. This example demonstrates that the cross-correlation 
technique works and is able to recover extremely faint emis- 
sion lines. 

4. LINE IDENTIFICATION 
4A. Is UDFj-39546284 at z= 121 

If the galaxy is at z ~ 12, as is suggested based on its strong 
photometric break between the JHuq and //i6o filters (Ellis 
et al. 2013; Bouwens et al. 2013), the line — if real — can be 
identified as Ly-a redshifted to z = 12.12. The rest-frame 
equivalent width of Ly-a would not be unreasonably high at 
> 170 A (accounting for the z = 12 Lyman break in //i6o); such 
values can be reached in young, low metallicity starbursts at 
high redshift (Schaerer 2003). However, such strong Ly-a 
emission might be unexpected early in the reionization epoch 
when the neutral fraction is high (Santos 2004; but see also 
Dijkstra et al. 201 1). More to the point, if the line is real it can 
account for most or all of the Hi^q flux, meaning that the pho- 
tometric break does not necessarily reflect a strong continuum 
break and the photometric redshift is unreliable. Therefore the 
line could also be a longer- wavelength feature at much lower 
redshift, a possibility discussed by Ellis et al. (2013). 



4.2. Low redshift solutions — [ O III] at z = 2.2 

Although we cannot exclude other low redshift identifica- 
tions such as [O II] at z = 3.28 or Ha at z = 1 .44, a likely possi- 
bility is that the line is [OIII]A4959,5007 atz = 2.19±0.01. A 
population of extreme [O III] emission line galaxies at z ~ 1.7 
was recently identified by van der Wei et al. (2011) from 
their significant line contribution to the /125 photometry in 
the CANDELS survey. Additional galaxies with [O III] rest- 
frame equivalent widths reaching 2000 A have been identi- 
fied in WFC3 grism spectroscopy by Atek et al. (2011) and 
Brammer et al. (2012b). The strongest emitters in van der 
Wei et al. (2011) have Jns-Hm - -1, with [O III] in the 
/i25 band. While such colors can mimic those of high-redshift 
dropout galaxies (Atek et al. 201 1), even more extreme equiv- 
alent widths at > 3000 A (rest- frame) are required to reach 
Ji25/i40-Hm > 1.4 observed for UDFj-39546284. 

4.3. Discovery of a Possible Analog: An Extreme [ O III] 
Emitter at z=l .605 

We have discovered a very bright example of such an ex- 
treme emission line galaxy in G141 grism observations of 
the EGS field (GO-12547, PI: Cooper). This remarkable 
object, EGS-XEW-1, located at a=14:17:58.2, (^=+52:31:35 
and shown in Fig. 4, is an "[O III] blob" at z = 1.605 with 
/Hi 40 = 21.5 and a complex morphology comprised of dif- 
fuse (A) and compact (B) components separated by V\ The 
exposure times and dithered pixel sampling for the Cooper 
et al. grism program are nearly identical to those of 3D-HST, 
and the spectra were processed with the 3D-HST interlac- 
ing software as described above. The diffuse component, 
which is extended over more than V' (8.5 kpc), has a com- 
bined observed-frame equivalent width of nearly 9000 A for 
the blended H/3 and [O III]A4959,5007 emission lines. The 
more compact component has EW=3000 ± 160 A. 

The broad-band SED of EGS-XEW-1 is compared to that 
of UDFj-39546284 in the right panel of Fig. 4. The spec- 
tral break that results from a line such as the one observed 
in EGS-XEW-1 is just large enough to satisfy the 2a limit on 
the JHuo-Hieo dropout color of UDFj-39546284. However, 
if UDFj-39546284 had exactly the same spectrum as EGS- 
XEW-1, the [O II]A3727 line would probably be detected in 
the bluer WFC3 HUDF12 photometry. Interestingly, EGS- 
XEW-1 has much redder UV colors than fx ~ observed 
for the high equivalent width starbursts at similar redshifts 
(van der Wei et al. 201 1 ; Brammer et al. 2012b), which, if also 
true for UDFj-39546284, could help explain the lack of z = 2 
Ly-a observed for that source (Ellis et al. 2013). While not 
a perfect match, EGS-XEW-1 provides a directly-observed 
route to a plausible low-redshift interpretation for the UDFj- 
39546284 photometry. 

5. DISCUSSION AND CONCLUSIONS 

We have described a deep HST grism spectrum of the can- 
didate z = 12 galaxy UDFj-39546284. Using the known emis- 
sion line morphology to increase the line sensitivity of the 
slitless grism spectrum, we detect a tentative emission line 
at A = 1.599/im with flux 3.5 ± 1.3 x 10"^^ erg s"^ cm'^ and 
observed-frame equivalent width > 7300 A. 

While the current observations do not conclusively forbid 
the z = 12 interpretation of UDFj-39546284, a number of inde- 
pendent factors conspire to suggest that the 2.7cr line is in fact 
real and that the true redshift of UDFj-39546284 is z 2.2. 
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Figure 4. WFC3/G141 spectrum and broad-band SED of EGS-XEW-1, an "[O III] blob" at z = 1.605 with extreme [O III] emission lines and colors similar 
to UDFj-39546284. The asymmetric line profile of component A is consistent with lines of H/5 and the [O III]A4959,5007 doublet as indicated; the lines are 
unambiguously resolved for component B, at the same redshift. The right panel shows the SED of the blob normalized to that of the UDFj-39546284 Hi^q 
magnitude. Non-detections for UDFj-39546284 are shown with 2a upper limits taken from Bouwens et al. (2013). Along with the G141 spectrum, ground-based 
photometry of the blob from the NEWIRM Medium Band Survey (NMBS) (Whitaker et al. 2011) is shown in blue, redshifted such that the [O III] line would 
have been observed at 1 .6 /im and scaled to account for the different NMBS-J3, JH^o, and //i6o passbands. The inset shows a color thumbnail of the blob created 
from HST ACS/y606, ACSZ/gn, and WFC3/JHuo images. 



First, the line is observed at just the right wavelength and flux 
to satisfy both the HUDF Hi^q detection and the >1.4 mag 
dropout at A < 1 .595 /im in the bluer HST bands. If the feature 
were observed just 4 pixels (93 A) bluer it would violate the 
constraints on JHuq by a magnitude. Second, Bouwens et al. 
(2013) demonstrate that without a significant contribution of 
Ly-a to the //i6o flux, the implied rest-frame UV luminosity 
for z ~ 12 is some 20 x higher than would be expected from 
the evolution of the luminosity function at z = 8-10. Third, re- 
cent deep spectroscopic surveys of photometric dropouts find 
a decreasing fraction Ly-a emitters with increasing redshift 
when the universe was increasingly neutral (Pentericci et al. 
2011; Caruana et al. 2012). Fourth, we show that spatially- 
extended objects with EW > 7300 A exist at z ~ 2. Finally, 
we find at least two [O III] emitters within 1000 km s~^ of the 
z = 2.19 solution for UDFj-39546284 in the 3D-HST cover- 
age of the full GOODS -South field, the nearest separated by 
69^' (570 kpc); the fact that other galaxies exist at this exact 
redshift increases the probability that the line is real and is 
[O III]A4959,5007. 

If the physical properties of UDFj-39546284 are similar 
to either the [O III] blob EGS-XEW-1 or to the high equiva- 
lent width starbursting dwarf galaxies studied by van der Wei 
et al. (2011), UDFj-39546284 at z = 2.2 would represent a 
new class of object 750 times fainter than the former or ~30 
times fainter than the latter. Scaling from the typical stel- 
lar masses of the van der Wei et al. (2011) galaxies, UDFj- 
39546284 would have a stellar mass of order 10^ Mq, similar 
to the mass of a single massive star cluster (e.g., Reines et al. 
2008). The extended //i6o morphology would then indicate 
the distribution of ionized gas surrounding the cluster. 

The most urgent question is whether the line is real and 
associated with UDFj-39546284 as opposed to coming from 
a nearby contaminating object spectrum or simply being a 
clump of positive noise fluctuations. If the line is real and 
is [O III] at z = 2.2, it should be possible to detect the 4959 



and 5007 A lines of the doublet separately, either in (much) 
deeper grism data or in deep ground-based spectroscopy. Ha 
may also be visible in the ^-band at a similar flux as [O III] 
(Fig. 4). Regardless of the reality of the emission line, we 
have shown that the broad-band SED of UDFj-39546284 pro- 
vides a reasonable match to a low redshift SED, specifically 
the observed spectrum of EGS-XEW-1, an "[O III] blob" at 
z = 1.605, redshifted to z = 2.2. This study demonstrates that 
the G141 grism on WFC3 can provide clean spectra with well- 
understood noise properties in deep integrations up to at least 
40ks (17 orbits). 

We gratefully acknowledge funding support from the Marie 
Curie Actions of the European Commission (FP7-C0FUND) 
and STScI grant GO-12177. 
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